A polar method for obtaining wave resonating quadruplets in finite depths

A polar method for obtaining wave resonating quadruplets {K₁, K₂, K₃, K₄} in the computation of nonlinear wave–wave interaction source term of the wave model is presented with results for both deep and finite water depths. The method first determines the end radial points of the locus equation for K₂, for each set of input wave vectors (K₁, K₃) on the symmetry. The locus of K₂ (and hence K₄) is then traced in the anti-clockwise direction starting with the maximum radial point on the line of symmetry. It is shown that when k₃>k₁, the number of points on the locus varies when the orientations of the input wave vectors are changed and reduces when the difference in the magnitude of the input wave vectors is increased. A significant advantage in this method is that the angular increment on the locus for K₂ can be kept constant.     

一种全自动提取居民地的方法

对遥感影像中居民地提取的各种方法作了一般性的阐述,并在吸收这些算法优势的同时,提出了一种全自动实现遥感影像中居民地目标的提取方法：先求取非最大值抑制后的梯度图,而后进行高斯模糊处理,以提高影像上居民地区域纹理特征一致性;通过增强大津法和形态学的方法全方位腐蚀与膨胀,完成面状地物居民地的提取;最后进行图象矢量化处理。论述方法的同时,展示和分析了实验结果。     

Singular points of the polarization tensor

The problem to compute the magnetic field above the chromosphere using data of the vector = B _t/Bt that gives the projected field direction can be solved with different approximations. The field of direction vectors is, however, not the only field accessible to observations. The Stokes parameters, which are components of the radiation tensor, can be measured at each point of the image plane. The directions of the eigenvectors of the radiation tensor define two mutually orthogonal systems of integral curves in the image plane. These families of curves have singular points, which are generally of different type than those of the vector field. When the morphology of H chromospheric fibrils are used to infer the topology of the magnetic field, a similar problem is met, suggesting that singular points should also be present there.     

Mineralogical and chemical evolution of the Ernest Henry Fe oxide–Cu–Au ore system, Cloncurry district, northwest Queensland, Australia

The Ernest Henry Cu–Au deposit was formed within a zoned, post-peak metamorphic hydrothermal system that overprinted metamorphosed dacite, andesite and diorite (ca 1740–1660 Ma). The Ernest Henry hydrothermal system was formed by two cycles of sodic and potassic alteration where biotite–magnetite alteration produced in the first cycle formed ca 1514±24 Ma, whereas paragenetically later Na–Ca veining formed ca 1529 +11/−8 Ma. These new U–Pb_titanite age dates support textural evidence for incursion of hydrothermal fluids after the metamorphic peak, and overlap with earlier estimates for the timing of Cu–Au mineralization (ca 1540–1500 Ma). A distal to proximal potassic alteration zone correlates with a large (up to 1.5 km) K–Fe–Mn–Ba enriched alteration zone that overprints earlier sodic alteration. Mass balance analysis indicates that K–Fe–Mn–Ba alteration—largely produced during pre-ore biotite- and magnetite-rich alteration—is associated with K–Rb–Cl–Ba–Fe–Mn and As enrichment and Na, Ca and Sr depletion. The aforementioned chemical exchange almost precisely counterbalances the mass changes associated with regional Na–Ca alteration. This initial transition from sodic to potassic alteration may have been formed during the evolution of a single fluid that evolved via alkali exchange during progressive fluid-rock interaction. Cu–Au ore, dominated by co-precipitated magnetite, minor specular hematite, and chalcopyrite as breccia matrix, forms a pipe-like body at the core of a proximal alteration zone dominated by K-feldspar alteration. Both the core and K-feldspar alteration overprint Na–Ca alteration and biotite–magnetite (K–Fe) alteration. Ore was associated with the concentration of a diverse range of elements (e.g. Cu, Au, Fe, Mo, U, Sb, W, Sn, Bi, Ag, F, REE, K, S, As, Co, Ba and Ca). Mineralization also involved the deposition of significant barite, K(–Ba)–feldspar, calcite, fluorite and complexly zoned pyrite. The complexly zoned pyrite and variable K–(Ba)–feldspar versus barite associations are interpreted to indicate fluctuating sulphur and/or barium supply. Together with the alteration zonation geochemistry and overprinting criteria, these data are interpreted to indicate that Cu–Au mineralization occurred as a result of fluid mixing during dilation and brecciation, in the location of the most intense initial potassic alteration. A link between early alteration (Na–Ca and K–Fe) and the later K-feldspathization and the Cu–Au ore is possible. However, the ore-related enrichments in particular elements (especially Ba, Mn, As, Mo, Ag, U, Sb and Bi) are so extreme compared with earlier alteration that another fluid, possibly magmatic in origin, contributed the diverse element suite geochemically independently of the earlier stages. Structural focussing of successive stages produced the distinctive alteration zoning, providing a basis both for exploration for similar deposits, and for an understanding of ore genesis.     

东亚地区春冬季与夏季大气环流异常相互关系的研究

利用1951~2000年共50年的北半球500 hPa月平均高度距平场资料和奇异值分解技术（SVD），重点对东亚地区季节间大气环流异常的相互关系进行了初步探讨。结果表明，东亚地区季节间大气环流异常存在着较为密切的关联，并且这种明显的非同步联系具有时空相关显著的特点，尤其是夏季大气环流异常与其前冬和前春大气环流异常的联系更为密切。当前冬和前春北半球东亚大槽和北美大槽及蒙古高压偏强（或偏弱），极涡偏弱（或偏强），中高纬度盛行经向环流（或纬向环流），以及低纬和热带地区高度正距平（或负距平）明显时，则对应于夏季东亚地区西太平洋副高和鄂霍次克海阻高强度偏强（或偏弱），位置偏南（或偏北），贝加尔湖阻高强度也偏强（或偏弱），但位置偏西（或偏东）的大尺度环流形势出现。当春季北半球大气环流具有上述特点以及夏季鄂霍次克海阻高和西太平洋副高强度偏强（或偏弱），位置偏南（或偏北），且极涡较弱（较强）时，则东亚地区秋季大气环流对应于蒙古高压加强（或较弱），西太平洋副高减弱（或加强），并向南和向东移动（或移动较慢），极涡向南扩散（或扩散减弱），大气环流向冬季过渡加快（或减慢）。另外，大气环流异常还具有一定的持续性特征。     

奇异值分解和奇异交叉谱分析方法在华北夏季降水诊断中的应用

利用1965~2000年华北5省市及相邻省73个地面观测站逐月平均降水场及北半球500 hPa高度场、北太平洋海温场资料, 采用奇异值分解 (SVD)、奇异交叉谱 (SCSA) 分析方法, 将华北夏季降水场分别与1月北半球500 hPa高度场、冬季北太平洋海温场进行了诊断分析, 得出奇异向量分布型及相互作用的耦合周期信号。在对前4对奇异向量的分析中发现, 华北夏季降水全区域为正距平时与1月北半球500 hPa高度场PNA遥相关型关系非常密切。ENSO对华北夏季降水的影响确实存在, 但华北夏季降水全区域为正距平时与冬季北太平洋ENSO关系并不明显。同时还找出了华北降水与北半球500 hPa高度、北太平洋海温场相互作用的关键区。在华北各型降水与高度场、海温场关键区相互作用的耦合周期中, 前者以准2~7年振荡为主; 后者则周期较长, 最短周期仍为准2年振荡, 最长周期为准10~11年振荡。以上结论为进一步研究华北夏季降水短期气候预测方法, 提供了参考依据。     

GRAPES全球奇异向量方法改进及试验分析

基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。      

统一方差分量估计公式

本文提出了一种适用于协方差阵奇异或非奇异、设计阵列满秩或降秩时的方差分量估计方法;其公式推导简单,形式统一,不需解线性方程组,同时可保证迭代计算方差分量的非负要求。作为特例还和Helmert法及MINQUE法作了比较。最后讨论了用真误差进行方差分量估计的计算公式,并给出了一个测距误差分析的实例。